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ABSTRACT 


A  technique  is  presented  for  Che  measurement  of  the  Grunelsen  param 
eter  (T)  of  a  solid  from  the  free  surface  velocity  Induced  by  pulsed 
energy  deposition.  The  technique  is  used  to  obtain  both  temperature- 
dependent  and  temperature-independent  Grunelsen  data.  The  measurements 
employ  a  velocity  interferometer  to  determine  the  free  surface  velocity 
of  solids  exposej  to  a  pulsed  electron  beam.  For  the  cases  where  r  is 
Independent  of  temperature,  the  r  values  are  obtained  from  the  maximum 
free  surface  velocity  as  a  function  of  incident  fluence.  For  the  temper 
ature-dependent  case,  T  is  extracted  by  an  analysis  based  on  the  differ¬ 
entiation  of  the  free  surface  velocity  with  respect  to  the  Initial 
temperature  of  the  solid.  The  analysis  tskes  into  account  the  effects 
of  both  the  finite  exposure  time  of  the  electron  pulses  and  the  delay 
time  of  the  interferometer.  Results  are  presented  for  Al,  Cu,  Ge  and  Si 
in  a  temperature  range  where  r  is  constant.  In  addition,  r  is  measured 
for  silicon  as  a  function  of  temperature  over  the  rerge  50°K  to  300°K, 
where  T  is  strongly  temperature-dependent.  In  all  cases  the  results  are 
in  excellent  agreement  with  the  thermodynamic  values. 
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I.  INTRODUCTION 


Whan  a  solid  la  exposed  to  pulsed  radiation,  either  photon  or 
particle,  the  resulting  rapid  heating  can  generate  thermoelastic  stresses 
which  then  propagate  through  the  medium.  The  magnitude  of  the  thermo- 
elastic  stress  developed  is  characterized  by  the  Grdnelsen  parameter  T, 
which  is  defined  in  terms  of  thermodynamic  parameters  of  the  solid  as 
T  “  B/pKfCy  where  6  Is  t he  volume  coefficient  of  thermal  expansion,  p  is 
the  density,  Kf  is  the  isothermal  compressibility,  and  Cy  is  the  specific 
heat  per  unit  mass  at  ccaatant  volume.  The  importance  of  T  in  describing 
the  thermoelastic  response  of  materials  arises  from  the  fact  that  T  is 
directly  related  to  the  rate  of  change  in  pressure  p  with  internal  energy 
density  E  at  constant  volume: 

r  ■  (H)v  (l> 

For  the  case  of  pulsed  energy  deposition,  the  target  sample  can  essentially 
be  considered  to  be  lnertlally  clamped  during  the  exposure  time.  Thus,  if 
T  is  constant,  the  change  in  pressure  £p(r)  at  £he  point  £  inside  the 
sample  due  to  an  energy  deposition  AE(r)  is  Ap(r)  -  T  AE(r) .  More  generally, 
if  T  is  a  function  of  temperature,  or  specific  energy,  the  rise  in  pressure 
is  given  by 

Eq+AE<?) 

Ap(r)  -  j  T(E)dE  (2) 

Eo 

where  Eq  is  the  initial  value  of  the  internal  energy  density. 

For  most  materials,  T  is  fairly  constant  at  room  temperature  and  above, 
and  their  thermoelastic  response  at  these  temperatures  is  adequately  char¬ 
acterized  by  a  temperature- independent  Griinelsen  parameter.  At  lower 
temperatures,  however,  the  variation  of  T  with  temperature  often  becomes 
significant.  In  some  materials  this  variation  is  quite  strong  with  F 
even  becoming  negative  at  sufficiently  low  temperature.  Therefore,  if 
one  wishes  to  describe  the  thermoelastic  response  of  these  materials  at 
low  temperatures,  it  is  necessary  that  the  temperature  dependence  of  T 
be  known.  In  this  paper  we  describe  and  demonstrate  a  technique  that  can 
be  used  to  measure  T  as  a  function  of  temperature.  The  measurement  is 
made  by  observing  the  free  surface  velocity  induced  by  exposure  to  a 
pulsed  electron  beam. 

It  has  been  shown  previously1 • 2> 3*1*  that  the  thermoelastic  response 
of  solids  can  be  utilized  as  a  means  of  measuring  T  for  solids  char¬ 
acterized  by  a  constant  f.  In  particular,  using  a  displacement  type 


*R.  B.  Oswald,  Jr.,  D.  R.  Schallhorn,  H.  A.  Eiscn,  and  F.  B.  McLean, 
Appl.  Phys.  Letters  1_3,  270  (1968). 

‘R.  B.  Oswald,  Jr.,  F.  B.  McLean,  D.  R.  Schallhorn,  and  L.  D.  Buxton, 
Appl.  Phys.  Letters  16^,  24  (1970). 

3R.  A.  Graham  and  R.  E.  Hutchison,  Appl.  Phys.  Letters  1_1,  69  (1967). 
‘‘R.  B.  Oswald,  Jr.,  F.  B.  McLean,  D.  R.  Schallhorn  and  L.  D.  Buxton, 
J.  Appl.  Phys.  42,  3463  (1971). 
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laser  Interferometer3**4  to  study  free  surface  motion  proved  to  be  a 
fruitful  and  direct  means  of  determining  I‘.  Displacement  Interferometry 
analysis  was  also  applied  to  the  study  of  the  thermomechanical  properties 
of  certain  materials  at  low  temperatures . 5  It  was  found  that  Introducing 
a  one-dimensional  thermoelastic  model  that  explicitly  accounts  for  the 
temperature  dependence  of  T  Is  indeed  necessary  for  the  description  of 
the  dynamic  response  of  these  materials.  In  this  analysis  the  known 
temperature  dependence  of  T  as  ascertained  by  thermal  measurements  of 
8,  Kf,  and  Cy  was  used  In  comparing  the  experimental  observations  with 
the  predictions  of  the  thermoelastic  model.  No  attempt  was  made  to 
measure  r(T)  with  the  displacement  interferometric  technique,  primarily 
because  the  reduction  of  displacement  data  to  values  for  T  requires  the 
unfolding  of  a  double  integral  and  the  resulting  accuracy  would  be 
relatively  poor. 

The  accuracy  of  the  results  for  Griinelsen  values  Is  improved  con¬ 
siderably  if  either  Induced  Ltress  or  free  surface  velocity  measure¬ 
ments  are  used,  for  then  only  a  single  integral  is  involved  in  the  re¬ 
duction  of  the  data  to  T  values.  Gauster6  recently  employed  a  stress- 
gauge  technique  to  relate  T(T)  with  measured  stresa.  In  the  present 
study  we  utilize  a  "velocity"  interferometer  of  the  type  described  by 
Barker7  to  measure  the  free  surface  velocity  induced  by  a  pulsed  electron 
beam,  and  we  describe  the  analysis  for  extraction  of  r(T).  We  note  at 
the  outset  that  the  precision  with  which  T  can  be  measured  by  inter¬ 
ferometric  methods  in  general  is  of  the  order  of  5%.  This  precision  is 
certainly  poorer  than  can  be  accomplished  by  careful  thermodynamic 
measurements  of  8,  KT,  and  Cv.  However,  interferometric  methods  are 
direct,  fast,  and  relatively  simple,  and  yield  results  within  the  pre¬ 
cision  necessary  for  most  engineering  applications.  Further,  by 
averaging  a  number  of  data  points,  the  precision  is  usually  improved 
considerably. 

In  the  next  section  we  describe  the  analysis  necessary  to  extract 
Griinelsen  values  from  free  surface  velocity  data.  In  essence,  the 
analysis  is  a  differential  one,  where  the  Gruneisen  parameter  is 
sampled  only  over  limited  regions  of  temperature,  or  internal  energy 
density.  Corrections  lor  the  finite  duration  of  the  electron  pulse  and 
for  the  delay  time  of  the  interferometer  are  included.  In  Section  III 
we  discuss  the  details  of  the  experimental  arrangement  and  procedures. 

In  the  final  section,  as  a  check  on  the  technique,  we  first  determine 
values  of  T  for  a  number  of  materials  over  an  energy  interval  where  T 
is  known  and  constant.  We  then  use  the  technique  to  determine  the 
temperature  dependence  of  T  by  measuring  T(T)  for  Si  over  a  temperature 
interval  from  50°K  to  300°K,  and  compare  our  results  with  the  results 
of  thermal  measurements. 


bF.  B.  McLean,  R.  B.  Oswald,  Jr.,  D.  R.  Schallhorn  and  L.  D.  Buxton, 
J.  Appl.  Phys.  4_2,  3474  (1971). 

*W.  B.  Gauster,  Phys.  Rev.  B  £,  1288  (1971)  . 

7L.  Barker,  Behavior  of  Dense  Media  Under  High  Dynamic  Pressure, 
Symposium  HDP  IUTAM,  Paris,  1967,  Gordon  and  Breach,  New  York,  N.  Y. 
p.  483  (1968). 
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II.  METHOD  OF  ANALYSIS 


The  measurement  of  the  GrUneisen  parameter  presented  here  Is  based 
upon  solutions  of  the  one-dimensional  thermoelastic  equation  of  motion 
which  Include  the  dependence  of  T  on  the  Internal  energy  densityB,  or 
temperature.  These  solutions  have  been  described  in  detail  previously;1**5 
thus,  we  briefly  review  only  the  results  which  pertain  to  the  present 
discussion.  We  consider  a  thin  sltb  oi  elastic  material,  initially  at 
temperature  T0  and  corresponding  specific  energy  E0,  exposed  to  a  pulsed 
electron  beam  of  uniform  flux  incident  along  the  x  axis.  Let  x  *  0  de¬ 
note  the  plane  on  which  the  bean,  is  incident  and  x  -  L  denote  the  plane 
of  the  rear  surface  where  velocity  measurements  are  made.  If  the  sample 
thickness  is  much  smaller  than  the  lateral  dimensions,  the  central  region 
of  the  sample  will  initially  respond  only  along  the  x  direction  in 
a  one-dimensional  manner5.  We  therefore  postulate  that  the  deposited 
energy  is  immediately  coupled  to  the  lattice  as  thermal  phonons  and  is 
characterized  by  an  x-dependent  deposition  function  AE(x)  (with  units 
energy  per  unit  volume).  Figure  1  shows  a  typical  energy  deposition  pro¬ 
file,  in  this  case  for  aluminum  6061  exposed  to  an  electron  spectrum  with 
a  mean  energy  of  1.96  MeV.  Note  that  the  energy  deposition  reaches  its 
peak  value  a  short  distance  xp  from  the  front  surface.  We  take  as  a 
convenient  measure  of  the  electron  beam  exposure  the  energy  deposition 
at  the  peak,  AEp  ■  AE(xp).  Also,  we  point  out  that  for  the  present 
experiments  the  sample  thickness  L  is  always  made  greater  than  the 
electron  range  in  the  material;  therefore,  AE(x)-K)  for  x  -  L.  We  assume 
that  the  temperature  variations  of  the  density,  p,  and  of  the  acoustic 
speed,  c,  are  small  in  comparison  with  the  temperature  variation  of  T 
and  thus  can  be  neglected.  In  addition,  we  assume  that  the  electron 
pulse  can  be  treated  as  a  square  wave  in  time  of  pulse  width  t0  and  we 
take  t  *  0  as  the  time  the  pulse  first  strikes  the  front  surface  of  the 
sample.  We  further  assume  that  attenuation  of  the  stress  pulse  during 
the  first  cycle  can  be  neglected  or  can  be  accounted  for  by  a  simple 
correction  factor  determined  empirically. 

In  this  paper  we  are  interested  in  the  motion  of  the  rear  surface, 
specifically  the  rear  surface  velocity,  which  from  the  results  of 
reference  5  is  given  as  a  function  of  time  by 

-L-ct+i0 

V  (L,  *■ )  -  (l/pc£.0)  I  dxK(x)  (3) 

L-ct 

waere  ■  ci0  and  K(x)  is  defined  by 
Ep(x) 

K(x)  -  J  T(E)dE  (4) 

Eo 


HW.  B.  Gauster, 
^R.  J.  Clifton, 


Phys.  Rev.  L87,  1035  (1969). 

J.  Appl.  Phys.  4^,  5335  (1970). 
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Figure  1.  Energy  deposition  profile  in  aluminum  6061  for  the 
1.96  MeV  electron  beam  exposure. 


for  the  range  0<  x  <L.  To  account  for  the  periodic  nature  of  the  solu¬ 
tion,  the  range  of  K(x)  is  extended  to  all  x  with  K(::)  being  defined 
for  other  values  of  x  by  the  periodic  conditions 

K(-x)  -  -K(x) 

K(x+2L)  -  K(x)  (5) 

K(0)  -  0 
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In  Eq.  (4)  Ep(x)  ■  E^  +  AE(x)  is  the  final  (x-dependent)  energy  density 
of  the  sample.  The  final  temperature  profile  Tp(x)  is  related  to  Ep(x) 
by  the  expression 


Tp(x) 

Ep(x)  -  p  1  CV(T)dT  (6) 

To 

Eq.  (3)  gives  the  calculated  velocity  history  of  the  rear  surface  is  long 
as  the  stresses  do  not  exceed  the  elastic  limit  of  the  material  and 
attenuation  is  negligible. 

Consider  now  what  one  is  actually  measuring  with  the  velocity  inter¬ 
ferometer.  In  particular,  we  must  account  for  the  effects  of  the  finite 
rise  time  of  the  interferometer.  In  figure  2  we  show  a  schematic10  of  the 
laser  velocity  interferometer  (after  ref  7).  We  follow  the  analysis 
given  by  Clifton.9  Essentially,  one  is  looking  at  the  Interference  be¬ 
tween  two  light  beams  One  beam  is  reflected  from  the  sample  rear 
surface  at  time  t  and  travels  directly  to  the  photomultiplier  tube  in 
time  t 1  at  the  photomultiplier  at  tine  t+t^.  The  other  beam  is  reflected 
at  time  t  -  Tq  and  is  subjected  to  an  optical  delay  of  duration  tq  thus 
also  arriving  at  the  photomultiplier  at  time  t  +  tj.  At  time  t  **  0  the 
sample  rear  surface  is  stationary  and  we  suppose  that  the  two  beat,  are 
adjusted  so  that  they  are  in  constructive  interference,  i.e.,  the  delay 
leg  is  adjusted  so  that  it  contains  an  integer  number  of  wavelengths  A. 

For  t  >  0  the  real  surface  begins  to  move,  and  the  phases  rf  the  two  light 
beams  are  shifted  relative  to  each  other,  the  phase  shift  being 

At  (t  +  t,)  -  C 2 / C L )  lU(t)  -  U(t  -  TD)I  (7) 

where  Cl  is  the  speed  of  light  and  U(t)  is  the  displacement  of  the  rear 
surface.  We  may  drop  the  time  tj,  as  it  merely  shifts  the  entire  time 
record  by  the  constant  value  tj.  Now,  every  time  the  phase  shift  is 
equal  to  an  integer  times  the  period  of  the  light  waves,  a  fringe  is 
detected  by  the  photomultiplier  tube-the  two  light  beams  are  again  in 
constructive  interference.  Thus,  the  fringe  count  N(t)  at  time  t  is 
given  by  N(t)  *  Ai(t)/iL  which,  by  (7)  and  rewriting  the  displacement 
as  an  Integral  over  velocity,  becomes 


N(t) 


(2/A) 


f  V(t')dt' 

t-TD 


(8) 


This  expression  shows  that  the  fringe  count  at  time  t  is  related  to  an 
"average"  of  the  velocity  over  the  time  interval  Tn,  centered  around  the 
time  t  -  td/2.  In  fact,  if  rp  is  sufficiently  small  we  can  apply  the 
mean  value  theorem  to  the  integral  in  (8)  to  obtain 

N(t)  =  2t0/A  V  (t  -  td/2)  (9) 


M.  J.  Berqer ,  in  Methods  in  Computational  Physics,  edited  by  B.  Alder, 
et  al  (Academic  Press,  New  York,  19f>3),  Vol.  1,  p.  135. 
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Figure  2.  Schematic  diagram  of  the  experimental  arrangement 
depicting  the  velocity  interferometer,  the  sample, 
the  calorimeter,  and  pulsed-electron  beam  exposure. 


L'ore  precisely,  however,  we  rewrite  (8)  as 

N(t)  -  2td/A  <V(t)>To  (10) 

where  <V(t)>Tl)  is  the  velocity  averaged  over  the  interval  tp  centered 
around  time  t,  and  is  defined  in  conjunction  with  eq.  3  by 

/•t+TD/2 

<V(t)>  -  1/ td  /  V(t’)dt' 

TD  t— td/2 


(ID 


l/pci0TD 


rt+lft/2 

t-T„/2 


dt' 


fL-ct'+i0 

/  dx'K(x') 

L-ct’ 


(It  is  understood  now  that  all  velocities  refer  to  the  rear  surface 
velocity.  Also,  we  have  again  shifted  the  time  base  by  the  constant 
value  td/2;  thus,  the  fringe  count  at  time  t+tj  +  Tp/2  is  related  to 
the  average  velocity  defined  by  Eq.  11  at  time  t.)  For  later  purposes 
it  is  convenient  to  define  auxiliary  variables  x  -  L-ct  and  x'  **  L-ct'; 
then  Eq.  (11)  becomes 


<V(t  ■  L-x/c)> 

tD 


rX+l d/2 

1/pc101d  J  dx' 

x-  Iq/2 


dx"K(x") 

(12) 
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where  In  this  expression  the  first  Integral  accounts  for  the 

velocity  averaging  in  the  interval  tj),  the  second  Integral  for  the  finite 
pulse-width  effect,  and  K(x)  accounts  for  the  energy  dependence  of  I'  via 
Eq.  (4).  In  essence,  a  record  of  the  fringe  count  as  a  function  of  time 
is  directly  related  to  the  velocity  history  of  the  reur  surface,  or  more 
precisely,  to  the  time  history  of  the  averaged  velocity  as  expressed  by 
Eqs.  (11)  or  (12). 


We  now  describe  how  we  may  extract  Gruneisen  values  from  rear  surface 
velocity  measurements.  Assume  for  the  present  that  1'  is  positive  for  all 
energies.  Then,  experimentally  the  method  reduces  to  measuring  the  maxi¬ 
mum  fringe  count  observed  during  the  first  pass  of  the  stress  wave  at  the 
rear  surface  as  a  function  of  the  Initial  energy  density  of  the  sample 
and  of  the  incident  electron  energy  dose.  That  is,  we  determine 


N  (E  ,  AE  )  ■  2i  /X  Max  {<V(t)>  (E  ,AE  )) 


max  o  p 


tD  0  P 


2t_/X  V  (E  ,AT  ) 
D  max  o  p 


(13) 


where  Vmax  (E0  ,AEp)  is  defined  as  the  maximum  averaged  velocity  observed 
for  the  Initial  specific  energy  E0  and  peak  energy  deposition  AEp.  If 
la  and  £[>  are  sufficiently  small,  then  we  see  from  Eqs.  (12)  and  (4)  that 
<V(t)>Tp  achieves  its  maximum  value  from  that  portion  of  the  stress  wave 

generated  by  the  beam  heating  in  the  vicinity  of  Xp,  the  position  of  the 
peak  energy  deposition.  Note  that  if  T(E)  is  positive  for  all  E>E0, 
then  the  quantity  K(x)  also  achieves  its  maximum  value  at  the  point  Xp. 
Letting  tm  be  the  time  that  maximizes  <V (t ) >tj)  and  using  Eqs.  (12)  and 
(4),  we  write  V^x  (E0,  AEp)  as 


(14) 

x+*o  Ep(x' ) 

/  dx’  /  T(E)dE 

x  E 

o 


V  (E  .  AE)  -  <V(t  )>  (E  .AE) 
max  op  m  Tq  o  p 


/ 


W2 


^  X 


vV2 


where  Xn,  -  L-ctm.  If  is  much  smaller  than  xp,  which  is  the  case  for 
the  present  experiments,  then  the  value  of  x^  depends  essentially  upon 
the  relative  sizes  of  Iq/2  and  xp.  If  J-d/2  <Xp,  then  for  all  practical 
purposes,  xm  coincides  with  the  peak  position  Xp.  On  the  other  hand, 
if  Id/2  >xp,  then  as  a  consequence  of  the  symmetry  conditions  (5),  the 
correct  choice  of  xm  is  easily  seen  to  be  tp/2.  That  is  to  say 

xm  =  Max  (xp,  fcp/2}  (15) 

where  the  approximately  equal  sign  is  used  because  in  principle  xm  could 
be  shifted  slightly  away  from  the  right  hand  side  of  (15)  due  to  smoothing 
effects  of  the  spatial  integrations  in  (14).  However,  for  the  work 
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reported  here  such  a  shift  is  entirely  negligible.  (There  could  con¬ 
ceivably  be  a  significant  shift  only  if  and  (or)  Iq/2  are  of  the 
order  of  or  greater  than  the  range  of  the  incident  electrons  in  the 
sample. ) 

There  are  several  ways  one  can  obtain  the  energy  dependence  of  T 
from  Eqs.  13  and  14.  One  direct  approach,  which  may  be  descr  bed  as  an 
integral  (or  non-local)  method,  is  simply  to  expand  T  in  a  power  series 
in  E,  retaining  as  many  terms  as  desired.  The  expansion  coefficients 
may  then  be  determined  by  fitting  to  the  measured  values  of  Vmax  (E0, 

AEp)  in  some  fashion.  Although  such  a  method  may  be  entirely  satisfactory, 
particularly  if  T  does  not  vary  too  rapidly  with  energy,  we  choose  to 
use  a  differential  (or  local)  method  of  analysis,  which  involves  expan¬ 
sions  of  T(E)  only  in  narrow  regions  of  energy.  This  method  offers  the 
advantage  that  the  errors  associated  with  the  expansions  of  F  are  re¬ 
duced  considerably  or,  equivalently,  only  a  few  terms  are  needed  in  the 
expansion  to  achieve  the  same  accuracy. 

- o  understand  best  the  motivation  for  our  differential  approach, 
note  that  in  the  limits  £o-*-0  and  (14)  reduces  to 


V  (E  ,  AE  )  *  1/ pc  J  T (E)dE  (16) 

max  o  p 

o 

where  we  define  Ef  -  Ep(xp)  »  E0  +  AEp.  (In  this  case,  of  course,  xm 
*  xp  and  AEp  *  AE(xp)).  Fir  this  limiting  case  T  (E)  can  be  expressed 
in  terms  of  (Eo,  AEp)  -  Vmax  (E0,  Ey)  in  either  of  two  obvious 

ways,  obtained  by  appropriate  differentiation  of  Eq.  16. 


/3V  \ 

(i) 

r(Ep) 

/  max  \ 

PC  \dEp  )  E0 

/  3V  \ 

(ii) 

r(Eo) 

I  max ) 

“  -PC  V  n  )  E_ 

o  r 

(17a) 


(17b) 


In  other  words,  we  can  obtain  T  as  a  function  of  energy  either  by 
measuring  VmaX  as  a  function  of  Ep  holding  the  initial  energy  constant 
or  by  measuring  Vmax  as  a  function  of  E0  holding  the  final  energy 
constant.  In  either  case  T(E)  is  determined  by  the  slope  of  the 
resulting  curve.  In  addition,  we  can  differentiate  (16)  with  respect 
to  E0  holding  AEp  constant,  obtaining  the  following  relationship 

(iii)  r(E0)  -  r(Ep)  -pc  AE  (17c) 

In  this  case  T(E0)  is  determined-provided  F(Ef)  is  known-by  measuring 
Vmax  as  a  function  of  Eq  holding  the  energy  dose  constant. 
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For  the  case  of  finite  t0  and  £[),  matters  are  not  quite  this  simple, 
of  course,  but  Eqs.  (17)  do  serve  to  Illustrate  the  general  idea  of  our 
procedure.  Before  discussing  the  general  case,  we  point  out  that  although 
in  principle  each  of  the  three  methods  is  valid  for  all  energies,  as  a 
matter  of  practical  convenience  they  lend  themselves  to  actual  analysis  of 
data  for  different  situations.  In  particular,  If  I'  Is  independent  of 
energy,  method  (1)  Is  the  most  straightforward  to  apply;  we  can  simply 
measure  Vmax  as  a  function  of  energy  dose.  (Note  that  in  this  case  (16) 
reduces  to  simply  Vmax  (Eo,  AEp)  ■  !’/pc  AEp).  Method  (i)  can  also  be  used 
to  determine  P(E)  at  high  energy.  Methods  (ii)  and  (ill)  are  better 
suited  for  extracting  GrQneisen  data  in  the  low  energy  region.  (Keep  in 
mind  that  in  practice  AEp  is  sizeable  and  can  be  much  greater  than  Eo). 
Method  (ill)  is,  of  course,  restricted  for  use  when  T  values  at  the  final 
energies  are  known,  whereas  method  (ii)  has  no  such  restriction.  However, 
if  P(Ep)  is  accurately  known,  method  (iii)  offers  greater  accuracy  in 
applications  to  be  discussed  further  below. 


The  procedure  we  employ  for  the  general  case  (non-zero  C0  and  £d)  Is 
as  follows.  First  we  write  down  the  expression  (14)  for  for  two 

sets  of  energy  parameters  (E0,  AEp)  and  (E0,  AF.p),  and  then  take  their 
difference 


V 

max 


V 


V 

max 


e’  x  +e./2  x+£ 

/o  r  m  D  .  r  o  ,  , 

r(E)dE  +  J  7*  /  T5-  (18) 

E  x  -t_/2  D  x  o 

o  m  D 


f(E)dE 


where  we  have  used  the  fact  that  the  first  energy  integral  is  independent 
of  x.  (We  are  setting  up  the  analysis  in  a  completely  general  way,  so 
that  our  expressions  may  be  applied  to  the  case  of  finite  increments  be¬ 
tween  the  energy  parameters  as  well  as  to  the  case  of  infinitesimal  incre¬ 
ments.  This  procedure  proves  convenient  in  application  to  actual  data. 

We  can  recover  the  differential  expressions  generalizing  Eqs.  (17)  by 
taking  appropriate  limits).  Eq.  18  is  exact  within  linear  thermoelastic 
theory.  To  proceed  further  we  write  the  limits  of  the  second  energy 
integral  as 


Ef(x’)  -  Ep  -  AEp  I 1-e (x 1 ) I 

E'(x')  =  E'  -  AE'  |l-e(x')| 
r  r  p 


(19) 


15 


where  we  recall  Ep  ■  Ep  (Xp).  Hence,  e(x)  la  the  normalized  depth  dose 
profile  auch  that  e(xp)  -  1,  l.e.,  e(x)  ■  AE(x)/AEp.  Next,  ye  expand 
r(E)  In  Taylor  sprlea  to  flrat  order  In  E  around  the  points  E0  ■ 

(E0  +  Eq)/2  and  Ep  “  (Ep  +  Ep)/2.  That  la.  In  the  flrat  Integral  we  let 

r<E)  ,  ((*,)  ♦  (If)^  (E-iF)  (20.) 

and  In  the  second  we  let 

C(E)  s  r(iF)  ♦  (ff)E.Jf  (E-Sf)  (20b) 


These  expansions  are  valid  as  long  as  P(E)  Is  approximately  linear  in 
the  energy  Intervals  E0  to  Eg  and  Ep  to  Ej-  and  provided  that  B0  and  ip 
are  small  compared  with  the  electron  range  in  th >  sample— so  that 
1  -  e(x')  <<  1  In  the  range  of  the  spatial  integrations  of  (18).  After 
some  straightforward  algebra  (18)  reduces  to 


V  (E  ,  AE  )  -  V  (E  ,  AE  ) 
max  o  p  max  o’  p 


where 


1/pc  r(£  )(e'-E  )  +  r(L)|E -e'-(ae  -ae’  )  la  I 

looo  rrr  PP 

-  r'(L)(AE +ae')/2  |(E  -e')  I,  -  (AE-AE')I2|1 
F  p  p  F  F  *  P  P  I 


r'(^F)  “  1  DE  1  E-i! 


(21) 


and 


/ 


X 

m  D 


vV2 


dx 

l 


x+B 


D  x 


/  °  I  l-e(x')l  n 


(22) 


The  factors  Ii  and  I2  include  the  corrections  due  to  both  the  finite 
exposure  time  and  the  finite  delay  time.  They  are  independent  of  the 
energy  parameters  E0,  E0,  AEp,  AEpj  they  depend  only  upon  the  values 
of  B0  and  Bq  and  on  the  shape  of  the  depth  dose  profile.  We  note  that 
if  B0  and  Bjj  are  much  less  than  the  electron  range,  then  I2  <<  I]  <<  1 
In  fact,  for  the  studies  reported  here  I;  is  entirely  negligible.  We 
retain  it  only  for  the  purpose  of  keeping  our  expressions  general. 
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Equation!  (21)  and  (22)  arc  our  basic  operating  equations.  The 
energy  dependence  of  the  Griinelaen  parameter  la  now  obtained  In  a  manner 
aiollar  to  the  methods  Indicated  by  Eqs.  17  for  the  simplified  case  of 
*o  ■  Id  “  We  consider  each  method  In  turn. 

(1)  Setting  E0  -  Eq  In  (21)  we  obtain  the  equivalent  of  Eq.  (17a) 


r(£p) 


pc 


V  (E  , 
max  o 


AE  )  -  V  (E  , 
p  max  o 

(AE  -AE')<1-I,) 


AE  ) 

_ E 


(23) 


Taking  the  limit  Ep  -  Ep  ■  AEp  -  AEp  -*  0,  we  have 


(24) 


As  pointed  out  earlier  this  method  can  be  used  for  determining  P  at 
higher  energies.  It  is  particularly  useful  when  T  Is  constant  for 
E  >  E0,  in  which  case  (24)  becomes 


r  - 

(ID 

r(£  ) 

o 


V  (E  ,  AE  ) 
pc  max  o  p 

1-1 i  AE 

P 

Pu.cing  Ep  -  Ep  in  (21)  yields 


-pc 


V  (E  ,  E_)  -  V  (E  ,  E_) 
max  o  F  max  o  F 


E  -E' 
o  o 


/AE  +AE*N 
+  r(Ep) ij  -  r'  (Ep)  [  f2-  pJ  i 2 


(25) 


(26) 


which  in  the  limit  Eq  -*  E0  corresponds  to  (17b) 


r(Eo) 


r(EF)ii  -  r’(EF)AEpi2 


(27) 


Hence,  if  T(E)  is  known  in  the  vicinity  of  Ep  we  can  obtain  r  values  at 
low  energy  via  Eqs.  (26)  or  (27).  For  the  case  of  finite  increments, 
however,  this  method  may  be  cast  into  an  alternative  form  which  elimi¬ 
nates  the  need  to  know  T(Ef).  If  instead  of  setting  Ep  ■  Ep  in  21),  we 
choose  Ep  =  Ep  -  (AEp  -  AEp)  I  \ ,  or  equivalently,  Eq  “  E0  +  (AEp  -  AEp) 
I I ,  we  obtain 
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-pc 


r<*0> 


V  (E  ,  AE  )  -  V  (E* ,  AE  ) 
max  o  p  max  o  p 

E  -E* 
o  o 


E 

o 


+  (AE  -AE')Ii 
P  P 


+  r' 


+  AE 


1  "  II  J 


(28) 


If  i0  and  IQ  are  sufficiently  small  compared  with  the  electron  range,  as 
is  the  case  for  our  studies,  the  last  term  in  (28)  is  small  compared 
with  the  precision  of  our  results  and  may  be  neglected. 


(iii)  The  expression  corresponding  to  (17c)  <s  found  by  setting 
AEp  ■  AEp  in  (21).  Noting  that  E0  -  Eq  ■  Ep  -  Ep  we  have 


r(£o)  -  r(feF)  -  r'(feF)  ae  u 


-  pc 


V  (E  ,  AE  )  -  V  (E  ,  AE  ) 
max  o  p  max  o  p 

E  -E* 
o  o 


I 

Passing  to  the  limit  E0  -*  E0,  this  becomes 

/  3V  \ 

P(E0)  ■  r(Ep)  -  r'  (Ep)  AEpI  j  -  pc 

Again,  r(£p)  must  be  known  to  evaluate  F(£0). 


,  AE 

o  '  p 


(29) 


(30) 


We  note  that  experimentally  it  is  convenient  to  present  the  data 
for  the  case  where  T  is  a  function  of  energy  in  terms  of  a  family  of 
plots  of  VmflY  vs  E0  for  a  series  of  AEp  values.  If  r(Ep)  is  unknown— 
or  we  prefer  not  to  use  it  if  known-then  Eq.  28  with  the  second  term 
neglected  can  be  used  to  obtain  r(E0).  To  minimize  the  error  asso¬ 
ciated  with  the  expansion  of  T  around  the  point  E0,  adjacent  curves 
(adjacent  AEp  values)  should  be  used.  On  the  other  hand,  if  r(Ep)  is 
accurately  known,  for  example,  by  either  method  (i)  cr  (ii),  then  the 
differential  expression  (30)  of  method  (iii)  can  be  applied  to  each 
plot  individually,  thereby  eliminating  the  expansion  around  E0.  (Note: 
the  differential  expression  (27)  cannot  be  applied  to  cn  individual 
plot  because  in  using  (27)  Ep  must  be  maintained  constant.  We  still 
have  an  error  due  to  the  expansion  of  T  around  Ef,  but  in  practice  r 
generally  varies  much  more  slowly  at  high  energy  than  at  low  energy, 
and  hence  the  errors  associated  with  an  expansion  at  high  energy  are 
smaller  than  those  associated  with  an  expansion  of  T  at  low  energy. 

In  addition,  when  l0  and  are  small  the  range  of  the  expansion 
around  Ep  is  correspondingly  small,  further  reducing  the  error  in 
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method  (111).  In  principle,  therefore,  method  (iii)  should  provide  the 
most  accuracy  for  temperature  measurements,  though,  of  course,  Eq.  28  of 
method  (11)  is  certainly  the  most  advantageous  for  use  when  one  wishes  to 
avoid  the  necessity  of  knowing  F  at  higher  energy.  We  shall  apply  both 
methods  to  reduce  the  data  in  the  piesent  study. 

Up  to  this  point  we  have  been  assuming  in  our  discussion  that  F  is 
positive  for  all  values  of  energy  under  study.  In  this  case  the  maximum 
rear  surface  velocity  results  from  the  stress  generated  in  the  vicinity 
of  the  peak  in  the  energy  deposition  profile;  the  time  tro  at  which  Vnjax 
is  achieved  corresponds  to  the  time  required  for  the  stress  wave  to 
propagate  from  the  position  xm  as  defined  by  Eq.  IS  to  the  rear  surface, 
l.e.,  tm  -  (L  -  xm)/c.  (See  discussion  preceding  Eq.  15.)  Now  we  inquire 
into  the  situation  for  which  T  may  be  negative  for  certain  energy  regions 
under  study.  This  situation  is  known  to  be  the  case  for  certain  materials, 
for  example,  the  covalent-bonded  cubic  crystals  Si,  Ge,  and  InSb,  for 
which  F  is  negative  for  temperature  T  *  0.2  0d  where  Bp  is  the  Debye 
temperature  of  the  solid.  First,  we  note  that  if  T  is  in  fact  negative 
everywhere  in  the  energy  region  E0  <  E  <  Ep,  then  our  analysis  goes 
through  exactly  as  before,  except  that  the  signs  associated  with  the 
rear  surface  motion  are  reversed.  The  displacement  of  the  rear  surface 
is  initially  in  the  negative  x  direction,  and  hence  V^x  corresponds  to 
the  maximum  velocity  achieved  in  the  negative  direction  during  the  first 
cycle.  Again,  V^x  results  from  the  stress  generated  in  the  vicinity 
x  'v  xp. 

If  T  changes  sign  in  the  interval  E0  <  E  <  Ep,  the  velocity  in 
general  would  not  attain  its  maximum  magnitude  at  the  time  tm*  Instead, 
the  time  at  which  the  velocity  is  maximum,  as  well  as  its  sign,  depends 
upon  the  explicit  nature  of  the  energy  dependence  of  F  and  the  energy 
deposition  function,  and  upon  the  values  ot  Ej  and  AE„ .  To  avoid  causing 
the  analysis  to  become  unduly  complicated  in  determining  the  values  of  T 
for  this  case,  we  can  simply  use  the  fringe  count  at  the  time  tm  at  which 
the  maximum  average  velocity  would  occur  if  F  were  everywhere  positive  (or 
everywhere  negative).  That  is,  in  our  analysis  we  use  as  defined 

by  Eq.  14  where  x„,  is  still  chosen  by  (15);  however,  we  realize  that  Vmax 
so  defined  may  not  correspond  to  the  actual  maximum  rear  surface  velocity. 
Of  course,  we  must  also  take  care  to  ascertain  the  correct  sign  of  Vmax 
in  this  case.  We  note,  however,  that  in  practice  T  is  generally  negative 
only  over  a  relatively  small  energy  interval  at  low  energies-if  T  is 
negative  at  all.  (Recall  that  the  specific  heat  of  solids  at  low  temper¬ 
ature  is  proportional  to  T^.)  Further,  the  values  of  AEp  that  must  be 
used  in  actual  experiments  are  sufficiently  large  that  the  final  energy 
values  Ef  are  rather  high  in  the  energy  region  where  T  is  positive.  Thus, 
it  turns  out  that,  even  when  Eo  ^  0,  the  maximum  rear  surface  velocity 
does  in  fact  usually  occur  at  the  time  tm  and  is  positive. 
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III.  EXPERIMENTAL  PROCEDURE 


To  establish  this  Grtineisen  parameter  measurement  technique  a  series 
of  materials  having  known  T  values  were  investigated.  Then  the  values  of 
T  obtained  by  applying  the  analysis  described  In  the  previous  nection  to 
the  measured  maximum  rear  surface  velocities  produced  by  pulsed  electron 
beam  exposure  are  compared  to  their  thermodynamic  values.  For  the  temper¬ 
ature-independent  measurements,  the  materials  chosen  were  Si,  Ge,  Al,  and 
Cu  having  a  range  in  T  from  0.43  to  2.2.  The  temperature-dependent 
measurements  were  carried  jut  on  silicon  over  the  range  In  temperature 
from  50  to  300°K. 

Figure  2  provides  a  schematic  diagram  of  the  experimental  arrangement 
depleting  the  electron  beam,  calorimeter,  sample  and  interferometer.  The 
electrons  emitted  by  the  field  emission  tube  of  the  flash  x-ray  machine 
stream  In  a  partially  evacuated  drift  tube  and  emerge  through  an  aperture 
with  a-2  mil  titanium  transmission  window.  As  previously  described*4  this 
arrangement  produces  a  stable,  reproducible  electron  pulse  whose  full 
width  at  half  maximum  is  approximately  30  nsec  with  a  fluence  uniformity 
of  ±10%  over  the  face  of  the  exposed  samples.  The  machine-charging  voltage 
for  the  exposures  reported  here  was  either  4.1  MV  or  5.0  MV  and  produced 
electron  spectra  with  mean  energies  of  1.96  MeV  and  2.50  MeV,  respectively. 

The  peak  energy  deposition  in  the  exposed  sample  was  determined  using 
thin,  in-line  copper  calorimeters  in  the  following  manner.  For  measure¬ 
ments  of  rear-surface  velocity  at  room  temperature,  the  average  dose  in 
the  calorimeter  was  measured  by  the  response  of  a  thermocouple  welded  to 
the  calorimeter.  The  energy  deposition  profile  in  the  calorimeter  and 
target  material  were  computed  for  an  incident  fluence  of  one  electron/ 
cm2  using  the  electron  transport  code  ZEBRA10  and  a  measured  electron 
spectrum*  (also  S.  Grayblll  and  G.  Ames,  private  communication)  cor¬ 
responding  to  the  operating  voltage  of  the  machine.  The  maximum  energy 
deposition  in  the  target  was  then  determined  by  multiplying  the  observed 
calorimeter  response  by  the  predicted  ratio  of  the  peak  dose  in  the 
sample  to  the  average  dose  in  the  calorimeter.  For  the  low-temperature 
measurements  with  silicon,  the  calorimeter  response  at  room  temperature 
was  normalized  to  correspond  to  the  peak  dose  appropriate  to  the  observed 
peak  velocity.  Then,  for  exposures  intended  to  be  at  a  constant  fluence, 
the  calorimeter  response  at  lower  temperatures  was  used  to  determine  the 
constancy  of  the  incident  electron  fluence  while  the  initial  temperature 
was  varied.  The  initial  energy,  E0,  was  obtained  from  the  measurement 
of  the  initial  temperature  using  the  Debye  specific  heat  with  a  Debye 
temperature  of  650°K. 

The  target  sample,  calorimeter  and  optics  were  supported  on  a  vibra¬ 
tion  isolation  table  independent  of  the  electron  beam  machine.  All  the 
samples  used  were  optically  polished  with  flat  and  parallel  faces.  For 
the  room  temperature  measurements  the  samples  were  disks  that  had  a  diam¬ 
eter  of  5.08  cm  and  a  mass  thickness  approximately  1.1  g/om2.  These 
samples  were  mounted  by  supporting  the  sample  around  its  edge  with  a 
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neoprene  O-ring  which  applied  e  slight  pressure.  Under  this  arrangement 
the  exposed  diameter  to  thickness  ratio  was  greater  than  10  to  1  for  all 
the  samples  exposed  for  the  room  temperature  measurements.  For  the  cry¬ 
ogenic  exposure  of  silicon,  the  samples  were  rectangular  slabs  that  were 
glued  to  the  cold  finger  of  a  variable  temperature  liquid  helium  Dewar 
flask  as  previously  described5.  The  cold  finger  had  an  aperture  to  permit 
observation  of  the  rear  surface  motion  of  the  sample.  To  insure  that  the 
effects  of  the  constraints  caused  by  bonding  the  sample  to  the  cold  finger 
would  not  be  observed  during  the  first  half  cycle  of  the  Induced  motion, 
the  diameter  of  the  aperture  was  made  greater  than  twice  the  thickness  of 
the  sample.  Thus  for  both  the  temperature  Independent  and  temperature 
dependent  measurements  the  observed  response  could  be  considered  to  be 
one-dimensional  for  at  least  the  first  half  cycle  of  motion. 

The  rear-surface  velocity  of  the  exposed  samples  was  measured  using 
a  velocity  interferometer  similar  to  the  type  developed  by  Barker.7 
The  Interferometer  employed  a  collimated,  single-mode,  6328&,  He-Ne  laser 
as  the  light  source.  The  laser  beam  was  directed  through  a  lens  which 
focused  the  beam  onto  the  rear  surface  of  the  sample  where  it  was  re¬ 
flected  along  a  separate  return  path  and  recollimated  by  a  second  lens. 

The  rear-surface  velocity  of  the  target  was  determined  by  the  Doppler 
shift  in  frequency  in  the  reflected  laser  beam.  The  frequency  shift 
was  measured  with  a  Mach-Zender  Interferometer  composed  of  the  beam 
splitter  and  two  retroreflectors.  The  Doppler  shift  is  detected  through 
the  interference  fringes  of  the  recombined  beams  as  described  in  Section 
II.  The  fringes  were  detected  with  a  RCA  Type  7265  photomultiplier  tube, 
whose  output  was  amplified  using  a  wide  band  amplifier  having  a  fixed 
gain  of  20  dB  over  the  frequency  range  from  200  Hz  to  200  MHz.  This 
signal  was  recorded  on  an  oscilloscope.  The  combined  rise  time  of  the 
photomultiplier,  amplifier  and  oscilloscope  was  approximately  12  nsec. 

The  maximum  "averaged"  rear  surface  velocity  was  obtained  from  the  re¬ 
corded  interferometer  trace  by  determining  the  maximum  fringe  count 
(occurring  at  time  tm)  and  using  the  velocity  Interferometer  relation¬ 
ship  given  T[,  Eq.  10.  For  the  measurements  reported  here,  the  difference 
in  optical  path  lengths  varied  between  8  and  52  m  with  corresponding 
delay  times  ranging  from  26  nsec  to  170  nsec. 

IV .  RESULTS 


A.  Energy- Independent  Grunelsen  Measurements 

To  verify  the  experimental  technique,  the  Grflneisen  values  for  Si, 
Ge,  Cu  and  A1  (6061)  wore  determined  over  an  energy  interval  where  T  is 
constant  using  Eq.  25.  The  results  were  compared  with  thermodynamic 
values  reported  in  the  literature. 1 1 t12  For  these  measurements  the 
values  of  VMX  (AE„),  AEp,  and  the  acoustic  speed  c  were  determined 
experimentally;  widely  published  values  were  used  for  the  densities. 

The  correction  term  Ij  for  the  effect  of  finite  pulse  width  and  inter¬ 
ferometer  delay  time  was  computed  in  accordance  with  Eqs.  22  and  15 
using  the  measured  values  of  t0  (■  30  nsec)  and  ip  along  with  the 

11  Karl  A.  Gschneidner ,  Jr.,  in  Solid  State  Physics,  edited  by  F.  Seitz 
and  D.  Turnbull  (Academic  Press,  New  York,  1970),  Vol.  16,  p.  410. 

1?K.  Bruqqer  and  T.  C.  Fritz,  Phys.  Rev.  157'  525  (1967)* 
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computed  deposition  profile  based  upon  the  measured  electron  spectrum. 
(See  fig.  1  for  aluminum.)  Vmax  was  measured  as  a  function  of  the  peak 
dose  deposited  In  the  material;  these  peak  doses  ranged  from  a  minimum 
of  3  x  10s  ergs/cm3  to  a  maximum  of  12  -  45  x  10p  ergs/cm3  depending 
upon  the  material.  Typical  results  for  these  measurements  are  shown 
for  A1  in  figure  3.  By  fitting  a  straight  line  to  the  Vmax  vs  AEp  data, 
the  errors  In  both  the  dosimetry  and  velocity  measurements  are  reduced 
through  an  effective  averaging  process.  The  acoustic  speed  of  each 
material  was  dett -mined  simultaneously  with  the  measurement  of  Vmax  by 
observing  the  time  required  for  the  stress  pulse  to  propagate  twice 
through  the  sample  and  then  combining  this  "round  trip"  time  with  the 
measured  thickness  of  the  sample. 

Results  are  shown  in  Table  I.  In  all  cases,  the  observed  values 
of  F  are  within  approximately  5%  of  the  tabulated  thermal  values.  In 
fact,  except  for  A1  there  is  agreement  within  2%.  However,  the  samples 
used  in  the  present  experiments  were  not  pure  A1  but  were  made  from  the 
aluminum  alloy  6061  (1.0%  Mg,  0.6%  Si,  0.25%  Cu,  0.25%  Cr,  97.9%  Al). 

B.  Energy-Dependent  Grtinelsen  Measurements 

The  changes  in  the  free  surface  velocity  of  silicon  resulting  from 
the  variation  of  F  with  specific  energy,  or  temperature,  are  illustrated 
in  Figures  4a  and  4b.  These  figures  show  the  velocity  interferometer 
traces  for  a  Si  single-crystal  sample-oriented  with  the  1 100 1  axis 
parallel  to  the  incident  electron  flux-exposed  to  receive  a  peak  dose 
of  10.7  x  108  ergs/cm3  but  at  two  initial  specific  energies  1.73  x  108 
and  16.5  x  108  ergs/cm3.  At  an  initial  energy  of  16.5  x  108  ergs/cm3 
(To  ■  244°K) ,  T(E)  is  very  near  its  high  energy  limit  of  0.43;  thus  the 
electron-beam  exposure  produces  a  velocity  history  which  is  character¬ 
istic  of  the  constant  high  temperature  value  for  T.  As  indicated  in 
figure  4a  the  peak  velocity  corresponds  to  a  fringe  count  of  2.2  at  680 
nsec.  A  marked  decrease  in  velocity  is  produced  by  lowering  the  initial 
energy  to  1.73  x  10®  ergs/cm3  (T0  *=  108°K)  while  maintaining  the  same 
peak  dose.  As  shown  in  figure  4b,  the  peak  velocity  to  correspond  to 
a  fringe  counc  of  1.0.  The  decrease  is  the  direct  result  of  the  re¬ 
duction  in  T  at  lower  specific  energy. 

The  effects  produced  on  the  measured  velocity  history  by  varying 
the  initial  energy  is  shown  more  explicitly  in  figure  5.  In  these  plots 
we  show  the  velocity  histories  as  obtained  directly  from  the  inter¬ 
ferometer  traces,  again  for  the  100  orientation  and  for  a  constant  peak 
deposition  of  11.1  x  108  ergs/cm3.  Therefore,  they  are  the  "averaged" 
free  surface  velocities  as  described  by  Eq.  11;  they  include  the  effect 
of  the  long  risetime  of  the  interferometer  due  to  the  large  optical 
delay  time  (tq  *  170  nsec)  used  for  these  measurements.  The  solid  lines 
in  figure  5  are  drawn  to  fit  the  experimental  data.  Two  effects  are 
observed  as  the  initial  temperature  is  decreased:  (1)  a  decrease  in 
velocity  at  lower  temperatures,  and,  (2)  changes  in  the  velocity  profile, 
particularly  over  the  early  portion  of  the  history.  This  change  in 
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MAXIMUM  VELOCITY  (CM/SEC) 


Table  I.  Comparison  of  Measured  Values  of  T  with 
Tabulated  Thermodynamic  Values  of  T  from 
the  Literature. 


Material 

c  (105  cm/sec) 

xD(nsec) 

X1 

pine  as 

^thermal 

Si  mu 

9.36 

170 

.047 

.42 

.43 

Ge  1111] 

5.95 

42 

.014 

.75 

.76 

A1  (6061) 

6.67 

42 

.006 

2.07 

2.17 

Cu 

4.87 

26 

.015 

2.00 

2.00 

Figure  4a.  Velocity  Interferometer  traces  for  a  Si  single  crystal 
exposed  with  the  | 100 |  axis  oriented  parallel  to  the 
direction  of  the  incident  electron  beam.  Both  traces 
correspond  to  a  peak  energy  deposition  of  10.7  x  108 
ergs/cm3  but  at  initial  temperatures  of  (a)  T0  =  244°K 
and  (b)  T0  «  108°K.  The  horizontal  and  vertical  scales 
are  200  nsec/major  division  and  0.20V/major  division, 
respectively. 
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Figure  4b.  Velocity  interferometer  traces  for  a  SI  single  crystal 
exposed  with  the  1 1001  axis  oriented  parallel  to  the 
direction  of  the  incident  electron  beam.  Both  traces 
correspond  to  a  peak  energy  deposition  of  10.7  x  108 
ergs/cm3  but  at  initial  temperatures  of  (a)  TQ  =  244°K 
and  (b)  T0  *  108°K.  The  horizontal  and  vertical  scales 
are  200  nsec/major  division  and  0.20V/major  division, 
respectively. 


shape  results  from  the  fact  that  each  portion  of  the  velocity  history 
corresponds  to  an  integral  of  T(E)  over  different  energy  intervals,  as 
can  be  seen  by  referring  to  Eqs.  4  and  12.  Therefore,  unless  T(E)  is 
constant  over  the  energy  intervals  of  interest,  the  changes  in  the 
velocity  history  with  temperature  should  in  general  be  non-linear. 

The  velocity  curves  of  figure  5  show  evidence  of  such  non-linear  behavior 
for  Si  over  the  specific  energy  interval  from  1  x  108  to  20  x  108  ergs/ 
cm3.  The  dashed  curves  in  figure  5  are  the  corresponding  velocity 
histories  calculated  from  Eq.  12  using  the  energy  dependent  T(E)  as  de¬ 
termined  from  our  measurements.  These  curves  will  be  discussed  in 
further  detail  below. 

To  demonstrate  the  use  of  the  analysis  developed  in  Section  II  for 
determining  an  energy-dependent  T(E),  the  GrQneisen  parameter  of  Si  was 
measured  over  a  range  in  initial  energy  density  from  1  x  108  to  24  x  108 
ergs/cm3  and  then  compared  with  thermodynamic  data13*14.  For  these 
measure  ients  single  crystal  samples  with  (111)  orientation  were  used. 

1  *D.  F.  Gibbons,  Phys,  Rev.  JL12,  136  (1958). 

1  *R.  II.  Carr,  R.  D.  McCammon,  and  G.  II.  White,  Phil.  Mag  1_2,  157  (1965). 
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Figure  5.  Comparison  of  measured  and  calculated  velocity  histories 
corresponding  to  -he  response  of  a  single  crystal  of  Si 
oriented  with  the  (100)  axis  parallel  to  the  direction  of 
the  incident  electron  beam  and  at  three  different  initial 
temperatures  (a)  273°K,  (b)  126°K  and  (c)  94°K.  The 
solid  curves  are  drawn  to  fit  the  data  points;  the  dashed 
curves  are  computed  from  Eq.  14  using  our  measured  energy- 
dependent  r. 


First,  we  used  Eq.  28  of  method  (ii)  with  the  last  term  omitted  to  obtain 
r(E),  in  which  case  it  is  not  necessary  to  use  Griineisen  values  for  higher 
energy.  Then  these  results  were  used  in  conjunction  with  Eq.  30  (method 
ill))  to  obtain  a  more  reliable  estimate  of  T(E)  at  lower  energy  densities. 
The  values  used  for  (E0,  6Ep)  were  determined  in  the  following  manner. 

Maintaining  the  peak  dose  at  a  constant  value,  the  maximum  rear  surface 
velocity  was  measured  at  several  initial  specific  energy  values  covering 
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Che  range  from  1  x  10®  to  24  x  10®  ergs/cm®.  These  data  were  then  fitted 
with  a  smooth  curve  thereby  giving  V^x  (E0,  AEp)  vs  E0  for  the  fixed 
value  of  AEp.  Again,  the  fitting  of  a  smooth  curve  to  the  data  effectively 
averages  out  experimental  errors  associated  with  dosimetry  and  with 
reading  the  maximum  fringe  counts.  The  procedure  was  repeated  for  six 
different  peak  doses  as  shown  by  the  family  of  plots  in  figure  6— those  in 
figure  6a  were  obtained  with  the  1.96  MeV  electron  beam  and  those  in  figure 
6b  with  the  2.5  MeV  beam.  This  family  of  curves  formed  the  data  base 
used  in  the  analysis.  For  the  analysis  p  and  c  were  assumed  to  be  constant 
and  equal  to  their  room-temperature  values.  Again,  the  value  of  Ij  was  de¬ 
termined  from  the  values  of  t0  and  td  (t0  “  30  nsec  and  tq  -  170  nsec)  and 
from  the  computed  deposition  profiles.  For  the  1.96  MeV  beam  I]  -  .047 
and  for  the  2.5  MeV  beam  lj  ■  .028.  (The  value  of  I2  is  of  the  order  of 
.001  in  both  cases;  hence,  the  second  term  of  Eq.  28  is  completely 
negligible. ) 


Figure  6a.  Maximum  velocity  vs.  initial  energy  for  silicon  single 
crystal  samples  with  I 111 I  orientation  exposed  to  six 
values  of  peak  deposition  AEp.  The  curves  in  (a)  were 
obtained  with  the  1.96  MeV  electron  beam  and  those  in 
(b)  with  the  2.5  MeV  beam. 
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Figure  6b.  Maximum  velocity  vs.  initial  energy  for  silicon  single 
crystal  samples  with  |111|  orientation  exposed  to  six 
values  of  peak  deposition  AEp.  The  curves  in  (a)  were 
obtained  with- the  1.96  MeV  electron  beam  and  those  in 
(b)  with  the  2.5  MeV  beam. 


In  figure  7  we  present  the  results  for  T  vs  E  using  Eq.  28.  To 
minimize  the  error  associated  with  the  expansion  of  T(E)  around  the 
point  E0  ■  1/2  (Eq  +  E0),  we  used  only  adjacent  curves  (for  the 

same  electron  spectrum)  in  obtaining  the  differences  Vmax  l(E0,  AEp) 

_  Vmax  (Eq,  AEp)|  .  More  explicitly,  for  given  (adjacent)  values  of 
AEp  and  AEl,  a  value  of  E0  is  chosen  and  the  corresponding  value  of 
vmax  (E0,  AEp)  is  read  from  the  experimental  curve  (fig.  6).  Then 
vmax  (Eo .  AEp)  is  found  from  the  AEp  curve  at  a  value  of  Eq  such 
that  the  requirement  Eq  ■  E0  +  (AEp  -  AEp)  Ij  is  fulfilled.  Then 
application  of  (28)  yields  r(E,j).  This  procedure  was  iterated  over 
the  entire  range  of  possible  E0.  In  figure  7  each  of  the  four  symbols 
represents  data  obtained  from  a  pair  of  adjacent  velocity  curves.  A 
visual-inspection  "best  fit"  curve  (solid  line)  drawn  through  the 
resulting  data  then  represents  our  experimentally  determined  ene  / 
dependent  T(E).  For  purposes  of  comparison  the  thermal  data  of 
Gibbons13  and  Carr,114  et  al  (dashed  curve)  is  also  shown  in  figure  7. 
As  can  be  seen  there  is  reasonably  good  agreement  between  our 
measured  T(E)  and  the  thermal  data  for  all  energy,  the  agreement 
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In  the  first  place,  i'(E)  varies  most  rapidly  with  E  ut  low  energy  and 
such  an  approx imat  Ion  is  expected  to  be  poorer  in  this  region  titan  at 
high  energy.  Secondly,  when  we  approximate  a  nonlinear  curve  with 
straight  line  segments  the  resulting  approximate  curve  will  lie  to  the 
concave  side  of  the  exact  curve. 

To  improve  our  measurement  of  1(E),  particularly  at  lower  energy, 
we  used  the  differential  expression  (30)  in  conjunction  with  the  above- 
results  from  method  (li)  for  T(E)  at  higher  energy.  In  tills  case  the 
analysis  was  applied  to  each  experimental  curve  in  figure  b  individually 
The  value  of  T  at  a  point  E0  was  determined  from  the  slope  of  the  maxi¬ 
mum  velocity  curve  at  E0  together  with  the  linear  fit  to  I'  at  Ep  «=  F.0  + 
AEp  in  accordance  with  (20b).  The  data  points  (up  to  six)  for  each 
value  of  E0  were  then  averaged  with  the  results  shown  In  figure  3.  Here 
the  error  bars  reflect  the  statistical  uncertainty  (not  experimental 
error)  associated  with  each  value  of  Ec.  Again,  we  compare  with  the 
thermal  data  (dashed  curve).  As  is  apparent  the  agreement  is  signifi¬ 
cantly  improved  at  lower  energy;  there  is  essentially  exact  agreement 
for  E  >  2  x  108  ergs/cm3.  The  large  discrepancies  at  the  very  low 
energies  are  probably  due  to  rather  large  errors  in  the  Initial  temper¬ 
ature  measurement.  One  point  of  note  is  t^at  we  are  quite  definitely 
measuring  negative  Criineisen  values  for  E  <  2.U  x  108  ergs/cm3.  Thus, 
our  technique  is  fully  capable  of  handling  this  region  even  though  the 
final  energy  densities  lie  well  up  into  the  region  where  f  is  positive. 

As  a  final  check  on  our  measurements  we  compute  the  velocity  time 
histories  corresponding  to  the  observed  ones  in  figure  5  using  Eq.  12 
and  the  measured  T(E)  curve  of  figure  8.  The  results  are  the  dashed 
curves  in  figure  5.  We  consider  the  agreement  entirely  satisfactory. 

It  demonstrates  the  validity  of  the  thermoelastic  model  which  forms 
the  basis  for  Eq.  12  as  well  as  the  reliability  of  our  technique  for 
determining  T(E).  We  note  that  the  computed  velocity  history  for  T0  = 
94°K  contains  an  initial  negative  velocity  region  corresponding  to  the 
fact  that  the  integral  over  T(E)  is  negative  in  the  tail  of  the  energy 
deposition  profile.  As  T0  is  lowered  even  further  this  negative 
velocity  region  increases  both  in  magnitude  and  in  extent  in  time. 

V.  DISCUSSION 


The  results  presented  in  this  paper  demonstrate  that  the  free 
surface  velocity  of  a  thermoelastic  solid  subjected  to  pulsed  energy 
deposition  can  be  used  to  determine  both  energy-dependent  as  well  as 
energy-independent  Grvinelsen  data.  The  experimental  results  further 
demonstrate  the  adequacy  of  the  model  of  thermoelastic  response  to 
account  for  the  effects  of  the  finite  pulse  width  of  the  electron  beam 
and  the  finite  rise  time  of  the  velocity  interferometer  used. 

This  technique  provides  a  complementary  procedure  to  those  using 
propagated  stress  or  free  surface  displacement  to  measure  T.  In  a 
previous  paper8  a  technique  was  outlined  (but  not  utilized)  to  measure 
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Figure  8.  Measured  energy  dependence  of  T  for  silicon  obtained  by 

applying  the  differential  expression  (30)  (method  iii)  to 
the  data  of  figure  6  and  using  the  measured  T(E)  values 
shown  in  figure  7  for  the  final  energy  densities  Ef. 

The  error  bars  reflect  the  statistical  uncertainty  of 
averaging  the  data.  The  thermodynamic  energy  dependence 
of  T  is  given  by  the  dashed  curve.13*14 


the  temperature  dependence  of  T  using  a  Michelson  interferometer.  How¬ 
ever,  the  free  surface  velocity  is  related  to  T(E)  through  a  single 
integral  over  energy  while  the  free  surface  displacement  is  related  to 
T(E)  through  a  double  integral  (a  spatial  integral  over  the  whole  of  the 
deposition  region  in  addition  to  the  energy  integral).  Therefore,  T(E) 
was  determined  from  free-surface  velocity  rather  than  displacement 
measurements. 1 3 
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As  transmitted  stress  Is  also  related  to  !'  via  a  single  Integral, 
the  use  of  stress  measurements  us  a  function  of  initial  temperature  is 
also  an  applicable  technique  for  measuring  [(E).  Recently,  Gauster6 
performed  such  measurements  in  aluminum  und  silicon  using  quurtz  gauge- 
pressure  transducers  bonded  to  the  samples.  The  quartz  gauge  measure¬ 
ments  offer  the  advantage  that  low  stress  levels  can  be  measured,  und 
hence  low  electron  fluences  can  be  used  with  correspondingly  small 
temperature  rises  In  the  samples.  However,  relatively  large-diameter 
gauges  must  be  used  to  achieve  sufficient  accuracy  ut  low  stress  levels, 
and  this  may  lead  to  difficulty  with  short  one-dimensional  read  times. 

On  the  other  hand,  with  the  velocity  interferometer  one  can  measure 
the  response  at  a  point  in  the  center  of  the  sample  rear  surface  and 
therefore  enjoy  the  full  one-dimensional  read  time  as  determined  by 
the  acousclc  speed  and  the  ratio  of  the  beam  area  to  the  sample  thick¬ 
ness.  For  many  high-fluence  electron  beam  facilities  currently  avail¬ 
able,  the  beam  area  is  typically  3  cm2.  Thus,  for  a  low  density 
material,  with  a  small  Grunelsen  parameter  and  low  acoustic  speed  the 
longer  read  time  attainable  with  the  velocity  Interferometer  may  be 
advantageous.  This  advantage  is  particularly  important  for  those 
cases  where  there  is  appreciable  attenuation  of  the  stress  pulses  as 
they  propagate  through  the  sample.  For  exposures  below  the  spall 
threshold,  it  is  possible  with  interferometric  techniques  to  observe 
the  motion  of  the  rear  surface  over  several  cycles  and  correct  for 
the  attenuation4.  An  additional  advantage  of  the  velocity  inter¬ 
ferometer  is  that  it  does  not  require  the  presence  of  a  second 
material  medium  to  measure  the  response.  For  high  temperature  studies, 
the  need  for  the  second  material  may  introduce  problems  associated 
with  the  bonding  of  the  quartz  gauges  to  the  samples. 

Additional  measurements  of  energy-independent  V  values  have  been 
carried  out  with  750  and  200  KeV  electron  beams  producing  peak  doses 
greater  than  1010  ergs/cm3,  thus  enabling  the  use  of  velocity  inter¬ 
ferometers  with  delay  times  of  6  nsec  and  less.  These  results  were 
equally  successful  in  demonstrating  the  utility  of  this  procedure. 

The  precision  of  an  individual  measurement  in  the  velocity  inter¬ 
ferometer  technique  is  limited  by  the  errors  in  dosimetry  and  measured 
velocity.  The  precision  of  the  dosimetry  technique  is  within  approxi¬ 
mately  ±10%,  while  the  fringe  count  can  be  read  to  within  approximately 
a  tenth  of  a  fringe.  Hence,  the  precision  of  the  velocity  measurement 
is  ±0.1  (X/2tq)  .  Thus  for  a  fixed  delay  time  the  relative  error  in¬ 
creases  as  the  velocity  decreases.  For  both  the  energy-independent 
and  energy-dependent  measurements  the  maximum  velocity  values  used  in 
the  analysis  were  determined  from  smooth  curves  fitted  to  the  data, 
giving  in  the  first  case  as  a  function  of  peak  dose  in  the  second 

case  as  a  function  of  initial  energy  density.  In  addition,  for  the 
energy-dependent  measurements  the  results  of  several  experimental  runs 
were  averaged  to  produce  T(E)  at  each  value  of  E.  In  this  manner  the 
random  errors  in  the  data  are  effectively  averaged  out  to  a  great  extent, 
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thereby  Inproving  the  accuracy  of  our  measured  values  of  I'.  As  we  have 
seen,  our  measured  values  of  T  all  lie  with  5%  of  the  thermodynamic 
values,  except  for  the  low  energy  region  In  the  temperature-dependent 
results  for  silicon,  where  a  relatively  large  constant  error  in  the 
initial  temperature  measurement  at  the  lower  temperatures  is  probably 
responsible  for  the  discrepancy. 
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